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We present results of phase shift for 1 = 2 S-wave hit system with the Wilson fermions in the quenched 
approximation. The finite size method proposed by Liischer is employed, and calculations are carried out at 
P = 5.9 (a -1 = 1.934(16) GeV from m p ) on 24 3 x 60 , 32 3 x 60 , and 48 3 x 60 lattices. 



1. Introduction 

Lattice calculation of phase shifts is an impor- 
tant step for understanding of strong interactions 
beyond the hadron mass spectrum. For scattering 
lengths which are the threshold values of phase 
shifts, there are already several studies in liter- 
ature. For the I = 2 irir scattering calculations 
have been carried out with the staggered , the 
standardfi||,[|] and the improved Wilson fermion 
actions Recent studies by JLQCD [Q] and by 
Liu et. al. || carried out the calculation at several 
lattice cutoffs and obtained the scattering length 
in the continuum limit. For the scattering phase 
shift, in contrast, there is only one calculation for 
7 = 2, which used lattice data to estimate an ef- 
fective 7T7r potential and a quantum mechanical 
analysis to calculate the phase shift from the po- 
tential 

We present a lattice QCD calculation of the 
1 = 2 S-w&ve wit phase shift using the finite-size 
method 0. We work in quenched lattice QCD 
employing the standard plaquette action for glu- 
ons at j3 = 5.9 and the Wilson fermion action 
for quarks. The number of configurations (lattice 
sizes) are 200 (24 3 x 60), 286 (32 3 x 60), and 56 



(48 3 x 60). Quark propagators are solved with 
the Dirichlet boundary condition imposed in the 
time direction and the periodic boundary condi- 
tion in the space directions. Quark masses are 
chosen to be the same as in the previous study of 
quenched hadron spectroscopy by CP-PACS ||] 
{i.e., m n /m p = 0.491,0.593,0.692,0.752 ) for all 
lattice sizes. Preliminary results of the present 
work was presented at Lattice'01 0]. 

2. Method 

The energy eigenvalue W p of an S'-wave tttt sys- 
tem with momentum p and — p in a finite periodic 
box of a size I? is shifted from twice the pion 
energy 2 ■ E p by finite-size effects. Liischer de- 
rived a relation between the energy shift AW P = 
W p — 2 • E p and the scattering phase shift S(p), 
which takes the form 



tan 5 (p) 
irLp 



-x p - A n x 2 p - B n x 3 p + Q{xp) , (1) 



*presented by N. Ishizuka 



where p 2 = n ■ (2n/L) 2 , (n = 0, 1, • • • 6), x p = 
AW P ■ {2E p L 2 )/(16ir 2 ) = 0(1/L), and v n , A n , 
and B n are geometrical constants. The scattering 
length is given by do = (tan 5(p)/p) p_>o- 

In order to obtain the energy eigenvalue W p 
we construct irit 4-point functions G p k(t) = 
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Figure 1. R p (t) and D p (t) for p 2 = 0. 

(0|r2 p (i)Slfe (0) |0). Here Sl p (i) is an interpolating 
field for the S*-wave 7T7t system at time t given 
by fl p (t) = 2 fl 7r(-fl(p),*M--R(p),*) wh ere i? is 
an element of the cubic group. In numerical cal- 
culations we construct the source operator f2fe(0) 
using the noisy source method with U(l) random 
numbers. 

Since the 4-point function G p k(t) contains 
many exponential terms as pointed out by Ma- 
iani and Testa jl(J , the extraction of energy eigen- 
values from G p k(t) is non-trivial. We solve this 
problem by the method proposed by Liischer and 
Wolf jy]]. In their method, one diagonalize the 
M pk (t ,t) = [G(i r 1/2 G(t)G(t )~ 1/2 ] P fc at each 
t where to is fixed at some small value. The eigen- 
values are given by X p (t) — exp(— W p (t — to))- 

In the actual diagonalization we have to cut 
off the set of momenta. Here we expect that the 
components of G p k(t) or M p k(to,t) for p,k < q 
are dominant for the eigenvalue X q (t) in the large 
t region, while the components p,k > q are less 
important. 

With this expectation, we calculate the energy 
eigenvalue for the ground state [n = 0) and the 
first excited state (n = 1) for V = 24 3 , and 
also that of the second excited state (n = 2) for 
V = 32 3 and V = 48 3 . The cut-off dependence is 
investigated by varying the number of momenta 
N > n for N = 0,1,2,3. 

3. Results 

In order to examine the effects of diag- 
onalization, we calculate two ratios defined 
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Figure 2. R p (t) and D p (t) for p 2 = {2-k/L) 2 . 

by R p (t) = G pp (t)/[G;(t)] 2 and D p (t) = 
K(-t)/[Gp(-to)/Gf(t)] 2 , where GJ(t) is the pion 
propagator with momentum p. If G pp (t) or X p (t) 
behaves as a single exponential function, we can 
obtain the energy shift AW P from the ratio R p (t) 
or D p (t) by a single exponential fit. 

In Fig. |l| we compare the ratios for p 2 = 0. 
The two-pion source is placed at t = 8. The mo- 
mentum cut-off is set at N = 1. The signals are 
very clear and diagonalizations do not affect the 
result. We checked the cut-off dependence by tak- 
ing N — 2 and confirmed that it is negligible. In 
previous calculations of scattering lengths 
the ratio Ro(t) was used to extract the energy 
shift A Wo- Our calculation demonstrates the re- 
liability of these calculations. 

We compare the ratios for p 2 — (2n/L) 2 in 
Fig. ||. The momentum cutoff is set at N = 1 
and N = 2. In contrast to the p 2 = case, the 
diagonalization is very effective. We can observe 
a convincing single exponential behavior only af- 
ter the diagonalization. The cut-off dependence 
is negligible. 

The analysis shown here, and additional ones 
for p 2 = (2ir/L) 2 ■ 2 for V = 32 3 and V = 48 3 , 
lead us to conclude that the momentum cut-off 
should be taken N > n for the energy shift AW P 
( p 2 = (2tt/L) 2 ■ n ). 

In Fig. H we plot our results for the scatter- 
ing amplitude T/(327r) = E p ■ (tan S{p)/p) at 
m^/nip = 0.491 and 0.593 obtained by substi- 
tuting our data for AW P into ([l]). In order to 
obtain the phase shift for various momenta at 
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Figure 3. Results for scattering amplitude. 

the physical pion mass, we extrapolate our data 
with the following fitting assumption : T j (327r) = 
A w ■ (ml) + A 2Q ■ (ml ) 2 + A 01 -{p 2 ) + A 02 -{p 2 ) 2 + 
An ■ {m 2 p 2 ). The fit curves are also plotted in 
Fig. |. 

In Fig. |^ we compare the calculated phase shift 
S(p) at physical pion mass with the fit above with 
experiments Our results for 5{p) are 30% 

smaller in magnitude than those of experiments, 
and our result for scattering length, flom T = 
-0.0272(16), differs from the ChPT prediction 
given by aom^ = -0.0444(10) |]|, but consis- 
tent with the result of JLQCD at the same lattice 
spacing, a m^ = -0.0300(31) (§. 

A possible cause of the discrepancy is finite lat- 
tice spacing effects. The JLQCD results for scat- 
tering length show sizable scaling violation [Q. 
Hence that of the scattering phase shift cannot 
be considered small. Further calculations nearer 
to the continuum limit or calculations with im- 
proved actions are desirable. 
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